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Valparaiso, Indiana 46383 Space Administration 
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Abstract 

This report describes the development of a software to aid in the design and analysis of AC 
power systems for large spacecraft. The algorithm is an improved version of Electric Power 
Research Institute’s (EPRI) harmonic power flow program, “HARMFLO”, used for the study of 
AC power quality. The new program is applicable to three-phase systems typified by terrestrial 
power systems, and single-phase systems characteristic of space power systems. The modified 
“HARMFLO” accommodates system operating frequencies ranging from terrestrial 60 hz to and 
beyond aerospace 20 kHz, and can handle both source and load-end harmonic distortions. 

Comparison of simulation and test results of a representative spacecraft power system shows 
a satisfactory correlation. Recommendations are made for the direction of future improvements to 
the software, to enhance its usefulness to power system designers and analysts. 

Introduction 

The advent of large spacecraft has resulted in their increased electrical energy consumption. 
Additionally, planned space exploration will require more capability for extended human and 
equipment support. These requirements are well served by the AC power system. However, the 
presence of switching devices in the converters at the source end of currently proposed, spacecraft 
AC system and their connected nonlinear loads can cause several problems in the power system. 
One of the more serious problems is the occurrence of harmonic resonance. This is a phenomenon 
whereby the harmonics generated by the switching devices and the nonlinear loads excite resonant 
modes caused by cancellation or near cancellation of inductive reactance by capacitive reactance. 

*1990 NASA/ASEE Summer Faculty Fellow at Lewis Research Center. 
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Among the problems caused by voltage and current resonance in a power system are insulation 
failure due to overvoltage, equipment malfunction resulting from high frequency and multiple zero- 
crossings of distorted waveforms, radio-frequency noise, and transmission and equipment current 
overloads. Generally, the generated harmonics tend to worsen the overall power quality. 

This report describes a harmonic power flow program, “HARMFLQ,” for identifying and 
evaluating resonance problems, among others, which can occur in either terrestrial or large 
spacecraft power systems. The program accepts both three-phase systems commonly found in 
terrestrial power systems, and single-phase systems characteristic of spacecraft power systems. 
The applicability of the program encompasses frequencies ranging from terrestrial 60 Hz to 
aerospace 20 kHz and beyond, and harmonic distortion at the source - and load-ends of the power 
system. 


In its present phase, the harmonic power flow program contains a mathematical model of a 
Mapham inverter [1], a key subsystem of a representative 20 kHz spacecraft power system used in 
developing and testing the algorithm. In consideration of criticality of payload weight, the salient 
features of lightweight and simple design, with a minimum number of components, make the 
Mapham inverter an attractive source of AC power for spacecraft applications. Inclusion of the 
inverter model in the algorithm can expedite required performance analysis of the spacecraft power 
systems during their design phase. The design and operation of the Mapham inverter is well 
documented in the literature [1-3]. Hence, only the particular aspects of the inverter operation 
which impacts its modelling are discussed in this report. 

Another essential subsystem of the representative power system used is the single-phase 
voltage controlled rectifier. As a potential source of harmonic current with attendant problems in 
the power system, the mathematical model for this rectifier is also developed. 

Simulated preliminary results obtained from the algorithm compare favorably with published 
test results. Recommendations included in the report point to future enhancements to the power 
flow program to make it more useful to power system designers and analysts. 
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The circuit diagram of a 


Mapham inverter is shown in 
Figure 1. The SCR's are 
switched at the desired output 
frequency, f a . This switching 
reverses the current through the 
capacitor, C r This current 
resonates at a frequency 
determined by the series 



Figure 1 - Circuit Diagram of a Mapham Inverter. 


combination of the inductor, L r , and C r . 
The reversing of the current at a time 
when the resonating current is in its 
negative half cycle will sustain the 
inverter’s operation as an AC source. 
Figures 2 through 4, which can be found 
in the publication by Sundberg, Brush, 



Button and Patterson [3], show the 
reversing currents, I L1 and 1^, though the 
inductors, L r /2, and the total resulting 
current through C r . It is important to note 
that unless the resonant frequency, f r , is 
the same as f s , a distortion occurs in the 
capacitor current waveform. This 



distortion results in harmonic currents 
which are injected into the power system 
by the inverter. Since the impedance of 
the load placed across the output 
terminals of the inverter will change f r , 
loading will affect the level of injected 
harmonic current. 



TM (IWX«) 


Figure 4 - Total Current Through Capacitor C r . 
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3 Modelling the Maoham Inverter 


To determine the injection currents that will occur for a certain power system 
configuration, it is necessary to find f r for the inverter under the desired loading. Under 
no-load conditions, the resonant frequency is given by 


fr 


1 

2n JLC 


When the inverter is connected to the power system, the system appears as an equivalent 
impedance of = R eq + JX eq connected across the output terminals of the inverter 
as shown in Figure 5. In the frequency domain, the input current resulting from a step 
input voltage has the form 


a m s m + ... + a, s + ^ 
b n s n + ... + b^s + b 0 


where, m < n and the coefficients of the two 
polynomials depend on whether the system 
equivalent is capacitive or inductive (i.e., 
whether = R eq - JX eq or Z eq = + JX eq ). The details concerning the derivation 

of these polynomials can be found in Appendix A. 

If the system equivalent is inductive then n = 5, and the coefficients of the 
denominator are: 



Figure 5 - Power System Equivalent. 
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i> 5 - ( Lj + L eq ) L r C r C s , 


■^4 = Roq^i C r C m , 


b 3 = (L 2 + L") C, + (l + * C.) , 


t >2 = ReqCs + 


RoqL t (C a + C z ) 


i\ = l + — * L<>< ? , and 


b = eq 

0 T 
■‘-‘m 


However, if the equivalent is capacitive then n = 6, and the coefficients of the 
denominator are: 

t>e = L 1 L I C eQ C I C g , 

•^5 — R oql- 1 1 C eqC I Cg ' 


■^4 - ^>1 C s + L z C z C g + L z C ai j C z + L z Cgq C e + 




b = r c c + ^iahi £gg 

^eq^eq’-s * } 
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A = C*q + + 


Ar ^ ^ w * At ^ 


2> x = , and 

An 


b °-± • 


Once the coefficients are computed, the denominator can then be factored into 


{s n + P n ) is + P n . x ) ... (S 2 + P 2 ) (s x + P x ) . 


The roots of this polynomial (i.e., the P's) represent the poles of the network. Since the 
network is primarily a series L-C circuit, at least two of the poles will form a complex pair 
with the form: 


( s z + a z + j a) z ) ( s 1 . 1 + a z - j a> z ) . 


Since n > m, this polynomial is said to be a proper rational polynomial and can be 
expressed in partial fraction form: 




+ ... + 




— + ...+ 


s r + o r -j(0 r s r + a x + jw r 


*i + Pi 
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The function can then be returned to the time domain by using inverse Laplace 
transforms. The inverse transform for the complex pair is 


cf-i 


KlQ + KL -8 
s + a-j'P s + a + j P 


2 Ke~ m * cos ( p t + 0 ) . 


(1) 


This yields a time domain response for the current through C r which is 


responses due 

I c ( t ) = J.e t sin ( o> r t + 0 r ) + to the other . 

V*r U A A. _ 

poles 



determined. Figure 6 shows one period of the current through C r . Also shown in Figure 
6 are the corresponding currents I L1 and I ^ which flow through each set of switched 
inductors. From this information, the capacitor current can be written as a function of 
time over five individual regions of the waveform. The resulting equations are: 
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I c (t) = -e “ 2 sin (» r t+w) 

for 0 s t < Ll ^ - l . 


( 2 ) 


J c ( t ) = e 




sin[to r t-7t(-^)] - 


e 2 sin(c o r t + ic) 


for lifLi s t < ^ , 


( 3 ) 


J C/ ( t ) = e 


-« tt- T *; T - t 




( 4 ) 


for-^ s , 


J c (t) = e 


-« ( t- - 


* 2 * ’ sin[o r t - n ( JL T7 i )] ~ e * <t T "* 2 * sin[o r 


( 5 ) 


for x t — ■—■ £ £ < , ancf 


J c ( t ) = - 

*-r 




sin[to r t-it(il^)] 
s t < x B . 


( 6 ) 


for -Lilli 
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Since the current is now mathematically represented as a continuous, repetitive time 
function, it can be represented by a Fourier series. The Fourier series of a repetitive time 
function is given by 

m 

f ( t) = c 0 + £ c D cos u<i)q t + d n sin no> 0 t 

fl -1 

where, 

ti-T 0 

C 0 = — / f U) dt , ( 7 ) 

V 


c n = — f f ( t ) cos (nti> 0 1 ) dt , ( 8 ) 

T ° t 

and 

d n = — f fit) sin (nti) 0 1 ) dt . ( 9 ) 


The Fourier coefficients, c n and d n , are found by integrating Equations (2) through (6) 
over one period of the waveform. It should be noted that the magnitude of the current 
is assumed to be 1.0 at this point. Once the coefficients are computed, they are then 
normalized by 
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Di L d'i = — L 8 j - nSi 


( 10 ) 


The coefficients are normalized because the magnitude of the current will not be known 
until a fundamental power flow is found. The required values of D 1 and 0 1 are then 
known. 

To apply these coefficients to the harmonic power flow algorithm, it is first 
necessary to understand the formulation of the harmonic power flow algorithm itself. The 
details of how this will be done are covered in Section 7 of this report. 
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4. Operation of the Single Phase. Voltage Controlled Rectifier 


The circuit diagram of a single phase, voltage controlled rectifier is shown in Figure 

7. Since the D.C. output voltage of the rectifier is the average of the "ripple" output 
voltage, V x , 
shown in Figure 

8, it can be 
controlled by 
changing the 
firing angle, a . 

This means that, 

by adjusting a, the average value of the "ripple" or the D.C. output voltage can be 
controlled within certain limits. The SCR's are triggered with a signal which is set by a. 
The value of a is determined by the difference between the desired D.C. output voltage 
and the actual D.C. output voltage. Once either of the SCR's is triggered, it will conduct 
until the voltage on the anode is less than the voltage on the cathode. This occurs when 
the input voltage, V in , drops below the output voltage, V out , atjs, the commutation angle. 



Figure 7 - Single Phase, Voltage Controlled Rectifier Circuit. 
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Figure 8 - Input Voltage to a Single Phase, Voltage Controlled Rectifier Circuit. 
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5. Modelling the Single Phase. Voltaae Controlled Rectifier 


When either of the SCR’s is triggered, the voltage which is applied to the input 
terminals of the rectifier is given by 

v^it) = v a sin (o,fc + o) u ( t) . 


This function becomes the following in the frequency domain: 


in 


(S) - V a 


r flflina «■ <i) g cos a 


The entire single phase rectifier circuit in the frequency domain is shown in Figure 9. The 
current source, CV(0 + ), is needed to account for any initial voltage that might be present 
on the capacitor when the SCR is fired at a. Since I in (s) is given by 


in 


is) 


Vj a (s) -V x (8) 

sL 


0 


V x (s) must be determined. This 
with the unknown voltage, V x (s). 


is done by writing Kirchhoffs current law at the node 



Figure 9 - Single Phase. Voltage Controlled Rectifier Circuit 
in the Frequency Domain. 
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( 12 ) 
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The result is 


sCV-ls) + 


V m ( s sin a + u> ,cos a ) ^ ^ 

V M (s) (s 2 + « 2 # ) 

R sL 


* CV[ 0*) 


K,(»)(sC* A.-i) • 


V,, ( s sin a ♦ « # cos a ) ♦ sLC ( s a * u a a ) V ( 0* ) 
bL(b 2 * a\) 


V x (s) 


££ [ V m ( 8 sin a + u^cos a) * sLC ( s 2 ♦ a t 2 ) V ( 0* ) ] 


( s ‘ ' ' Tc) ' 


If Equations (11) and (13) are substituted into Equation (12), the result is 


V m (ssina ♦ o.cosa) -^ [ V. ( s sin o + o.cos « ) + sLC ( s* + «*. ) V<0*>] 


( s 2 + ) 


I mis) = 




SL 


s|s + -gg |[ V a ( ssina + a> # cos a ) ] - s ( s 2 + o> 2 # ) ^(0*) 

SL ( sI * ~Kc * Tc ) ( s * + “ a «> 


| s » | [ s sin a + u, cos a ] - ( s 2 + &> 2 , ) V ( 0* ) 


( 


s 2 + + 


AC LC 


j ( s 2 + u> \ ) 


( 13 ) 


13 



This can now be expressed in partial fraction form: 


inis) = 


*1^1 


jqz-Oi 



f 1 \ z 

V jlc ' 

i jg] 


2RC 


& r + 


x ,j /¥W 

r* /nr ^ n 


2RC 


+ * 2^ 0 2 + ^^- Q 2 
S-JU, S+j«, 


where, 


*1 L e i = 


s + 


2 tfC 




\ RC/ 

Ti»(s) 

2 


L. ■ j ^ *C '( «c) 


2 J?C 


and 


jc 2 Z6 2 = [s- jo s ] J^fs) 1,.^ • 


To simplify the notation, let 


o 


l 

2 RC 


and A 
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Once the complex coefficients, K 1 and K 2 , are computed, the transformation back Into 
the time domain is found using Equation (1). Now that the input current is known in the 
time domain, Equations (7), (8), and (9) can be used to find the Fourier coefficients of the 
current. It is important to note that, because the total current is known at this point, 
normalization of the rectifier coefficients is not needed. 

These coefficients, like those found for the Mapham inverter, must be incorporated 
into the harmonic power flow program. To accomplish these goals, the following is a brief 
summary of the development of the harmonic power flow program. 
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6. The Harmonic Power Flow Algorithm 


The harmonic power flow algorithm used in this work is the one developed by 
Heydt and Xia [4,5] and later modified by Grady [6]. The algorithm uses a standard 
Newton-Raphson formulation for a power flow program. The basic algorithm is: 

1. Select an initial solution vector for all bus voltages, [V {0J ], where 

element j is | | and e i or the complex voltage at bus j. 

2. Using the voltage vector, calculate 

fa = p ^scheduled _ p ^calculated 


4 0, ■ Oj 


scheduled _ q calculated 


If A P and A Q for ja[l busses are < e , a specified tolerance, a solution 
has been reached; and, the iterative process is stopped; and, the 
output is printed. 

3. Compute the elements of the Jacobian matrix which are defined by 


iJi,] 


dp i 

dP, 

dVj 

v jMj 

BQi 

so, 

CD 

V 5 ; 
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4. Using [A P, A Q] T = [J ] [A V, A 6 ] T , calculate A V and A 6 . 

5. Update the solution voltage vector, 

[V {i} ] = [V {M} ] + [AV, AS] 7 . 

6. Repeat Steps 2 through 5. 

This basic power flow algorithm was modified by Heydt, Xia and Grady in the following 
manner. 

1. The basic algorithm described above is run on the system for the 
fundamental frequency only. 

2. The solution voltage vector obtained in Step 1 is used as the initial 
value in the harmonic power flow. Therefore, the initial voltage 
vector is 


y(l>‘ 
tr<a) 
V ( 3 ) 


y (a) 


where, V (i) represents the bus voltage for the i m harmonic and • 
are the firing and commutation angles for the converters. 
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3. Using the solution vector, calculate 


A Pj = 


d scheduled 


o calculated 


A Qj 


q ^ scheduled _ £^calculat*d 


^ j- ( i ) - j ( i ) scheduled _ y ( i ) caicuiated 


( 14 ) 


where, l (l) scheduled is the scheduled i th harmonic current at each bus, 
and |® calcula,ed j S the calculated i ,h harmonic current at each bus. 
It should be noted that the scheduled current for all busses except 
harmonic source (i.e., nonlinear devices) isO. If AP, AQ, and A I are 
< e , a specified tolerance, for Ml busses, the solution has been 
reached; and, the algorithm is finished. 

4. Calculate the elements of the modified Jacobian matrix. The 
Jacobian matrix is modified by formatting it as follows 


' J n) 

J< 2) 

. . . jin) 

0 

J*Q ( 2 , 1 ) 

rpQ ( 2 , 2 ) 

♦ • • j ->Q(2,h) 

H i2) 

(h, 1 ) 

fj*Q {h,2) 

. . . •j'Qihth) 

H ih) 

TO* 1 ' 1 * 

y (1,1) 

. . . y*£(l,h) 

H m 
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where, 


and 



dP; < m > 

t 

1 

... fin) _ 

J U 

d Vj m 
d Oj ( m ) 

Vj d bj (m) 
d Oj ( m ) 


dVj (m) 

Vj d bj < m > 

TGjj ( m ’ n > - 

dRe(l, ( m >) 
dVj w 

dim (I, ( m b 

JtpV— 

dRe (1/ < m b 
dbj (") 

dim Of ( m b 
d bj (n} 

Hjj ~ 

3 Re Oj (m b 
daj 

dlmOi (m) ) 

da] 

dRe Oj (m) ) 
dPj 

dlmOi (m) ) 
dP] 


5. Using [A P, A Q, A I] T = [J] [A V, A 8, A a, A p] T , calculate A V, A 8, A a, 

Ap. 


6. Update 

[VW.SW.aW.pW] 1 « [Vt M >, St'' 1 ), p{ M ^ T + [AV, A8, Aa,AP] T . 

7. Repeat Steps 3 through 6. 

The modifications introduced in the basic power flow make the HARMFLO different from 
others which use the fundamental power and harmonic current responses of the 
nonlinear devices to solve the voltage levels within a power system [7]. Furthermore, 
the HARMFLO differs from the Alternate Transients Program (ATP) version of the 
Electromagnetic Transient Program (EMTP) [8]. The ATP version is used in analysis in 
which the nonlinear loads are represented by harmonic current injections at desired 
nodes within the power system. 
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7. Modifying the Harmonic Power Flow Algorithm 


Since three phase, synchronous machines used in terrestrial power systems are 
designed to produce an undistorted 60 Hz sine wave under normal operating conditions, 
these electric power producing machines do not inject harmonic currents into the power 
grid. Because of this fact, the harmonic power flow program described above does not 
include the ability of the electric power sources to inject harmonic currents into the 
electrical system. Inverters, particularly the Mapham inverter, are, however, large 
producers of harmonic injection currents. These inverters inject the harmonic current at 
the bus where they are connected to the system. This change in sources required a 
modification to the harmonic power flow program. 

As previously stated, the Fourier coefficients of the capacitor current, l & , are 
computed by inserting Equations (2) through (6) into Equations (7) through (9) and 
integrating over the appropriate ranges. The resulting complex coefficients are then 
normalized by Equation (10). This process is summarized by 


IjU) = Re 


£ (£>/"> ze/ n >) e Jm, « c 


n* 1 


= Re e ja ° c + D/ 2, Z0/ 2 > e i2 ^ z + ...) 


Ij (t) noriMlized = Re ( 1 Z 0° e jUat + (£^)z (0 i < 2) -nQj' 1 ' ) e j2 * ot + ...) 


= Re (1Z0° e ju ° c + D i (2>/ Z0 2 / e J2 “ #t + ...) . 


Therefore, the magnitudes of the actual injection currents at every inverter bus for each 
harmonic can be determined by multiplying the magnitude of the fundamental source 
current at the bus by the appropriate Fourier coefficient. Thus, the scheduled injection 
current for an inverter bus j during the solution of the harmonic power flow can be 
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determined as follows by using the fundamental current magnitude and angle, lj (1) and 
iK j (1) , which is found in Step 3 of the modified algorithm described in Section 6: 

Ij ( t ) • ehedul9d * ( Ij (1) Zijr j (1) ) {Ij ( t) nornaliz0d ) 

= ( Jj (1> Zf/ 1 * ) [Re ejw 0 t + Z>/ 2>/ Z0/ 2)/ e J2 *' #t + - 

= Re [ Jj ( 1 > Zi|r j <l) e Jitot + J/ l> (D^') L (e/ 2,/ + ni|r J (l) ) e j 
= J?e [ (1) Z ^ (1> e J<1>ot + J/ 2 , Zt|r i < 2> e j2Wot + ...] . 

The injected harmonic currents are then placed in the appropriate locations in Equation 
(14). In order to place the harmonic injection currents into Equation (14), it is necessary 
to represent them in rectangular form (i.e., X + jY). Therefore, the current at each 
harmonic is found by 


j^(in) scheduled _ j (m) cos ( ^ .(®) ) 


= (j J . (1) ) (D j {m) ‘) cos(Q j (lB)/ - nQj ll) ) 

= (i j {1) ) ( ^Cj <in> ' 2 + ) cos tan' 1 ^ dj ( — j 


-ne/ 1 ’ 


and, in a similar manner, 


t cm) scheduled _ / j (i) ) ( - dJ jn) ' 2 ) sin tan' 1 / —!■ —r I - JiO x (1 

“‘‘J# Xm J V J J I c ^<*) ; ^ 
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The harmonic power flow program does model nonlinear devices such as rectifiers. 
The problem is that these rectifiers are three phase. Three phase rectifiers do not 
mathematically function like a single phase rectifier. For this reason, modifications in the 
harmonic power flow code are necessary. The modifications needed to successfully 
incorporate single phase rectifiers into the harmonic power flow code are simply to place 
the various harmonic injection currents described in Section 5 of this report into the 
scheduled current vector: 


T (m) scheduled 
-L j. Re 


and 


T (in) scheduled 
1 j.Im 


Once this is done for all Mapham inverters and single phase rectifiers, the 
harmonic power flow program is allowed to proceed per the algorithm described in 
Section 4. 
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8. Examples 


To demonstrate the practicality of the theories presented in this report, the 
following examples are offered. The first example will investigate the voltage distortion 
experienced in the power system depicted in Figure 10. This example is drawn from data 
obtained from the 
General Dynamics 
Corporation (GDC) 20 
kHz Testbed which was 
developed to study the 
proposed 20 kHz, A.C. 
power system. The 
inverter is a Mapham inverter set which will produce 440 at 20 kHz. This inverter will 

have the following electrical characteristics: 

L r = 16 \lH . 

C r = 1.71 F , 

C s - 2.0 [if, 

L 1 = 1.8 |Jl H , 

L m = 1.0 mH , and 
f s = 20 kHz . 


1 — 0 — 

Transmission Line (50 meters Ion*) — 


Inverter 
1 


A.C. Load 

— m — 


Figure 10 - General Dynamics Corporation (GDC) 20 kHz Testbed Power 
System. 


The transmission line is a 50 meter length of stripline designed cable which has the 
following electrical characteristics: 
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Resistance 


Inductance 
Capaci tance 



The load is standard A.C. load which will be varied from 2 kVA to 10 kVA. The load will 
maintain a constant power factor of 80% lagging throughout this example. Figure 1 1 
shows the data formatted for the harmonic power flow program. The load was varied 
from S = 0.16 + j 0.12 p.u. to S = 0.80 + ] 0.60 p.u. The resulting Total Harmonic 
Distortion (THD) of the bus voltage, defined as 



A ( (n) ) 2 

v (i> 

2 v i 


was used to 
generate the 
graph shown in 
Figure 12. Note 
that the results of 
tests conducted 
by Sundberg, et. 
al. are also 
presented on the 
graph in Figure 

12. Examination of the graph in Figure 12 reveals that the harmonic power flow does 
yield results which are very close to the values obtained experimentally. It should also 
be noted that the results obtained from the harmonic power flow program are 
exceptionally sensitive to the impedance values. It was found that very small changes in 
the data will cause relatively large changes in the output. 



Figure 11 - Input Data from the GDC 20 kHz Testbed for the Harmonic Power 


Flow. 
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The second example will also use the system described in Figure 10. This time the 
switching frequency will be varied to get the ratio f r /f, to vary from 0.6 to 0.95. The data 
for this example remains the same as that found in Figure 11 with the exception of the 
switching frequency. This value is varied as described above; and, the resulting THD of 
the bus voltage is compared to that measure by Sundberg, et. al. in Figure 13. 


p . £ . - 0.8, lagging 



KVA 

Invtrttr uata * 2.0 uf Sariea Capacitor 


Figure 12 - Comparison of Results of Bus Voltage THD versus KVA Loading from 
Harmonic Power Flow Program and Test Data from the GDC 20 kHz Testbed. 
Current THD versus KVA Loading is also Included. 
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Figure 13 - Comparison of Results of Bus Voltage THD versus f, n from Harmonic 
Power Flow Program and Test Data from the GDC 20 kHz Testbed. 
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9. Conclusions and Recommendations 


This report for the 1990 Summer Faculty Fellowship Program documents the 
development of software which will benefit the design and analysis of large A.C. 
spacecraft power systems. Early results indicate that the models developed are quite 
adequate for the intended goals. This is demonstrated by the acceptable correlation of 
the results exhibited in Section 8. Although the results are not exact, the discrepancies 
are relatively small. Given the small discrepancies and the limited test results from the 
GDC 20 kHz Testbed for comparison, further comparisons are necessary before accurate 
conclusions concerning program accuracy can be drawn. Once the GDC 20 kHz Testbed 
is made operational in its new location, an organized set of tests should be devised which 
will thoroughly evaluate the models developed in this report. This evaluation will help 
refine the models. 

As a result of the work cited in this report, the following areas are recommended 
for further investigation and development: 

• complete voltage control for voltage controlled rectifier, 

• model series and parallel operation of Mapham inverters, and 

• model bidirectional receivers. 

Implementation of the voltage controlled rectifier model discussed in Sections 4 and 
5 was started during this time frame. The model, however, was not completed due to 
time constraints. The model, which is presently incorporated into the harmonic power 
flow program, is a single phase, full wave rectifier. The model for this device is exactly 
that of the single phase, voltage controlled rectifier with the exception that the voltage is 
determined by the system configuration and not controlled to a desired level. The work 
required to implement voltage control appears to be straightforward. 
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Mapham inverters are operated in series to increase the total output voltage of the 
source. The modelling of this operating configuration must be investigated so that this 
mode of operation can be included in future design reviews. 

The bidirectional receiver unit is a device which is proposed for use on space 
vehicles. Modelling this device in the harmonic power flow program is important for a 
complete analysis of any proposed power system. For this reason, it is recommended 
that work be conducted in this area as well as those mentioned above. 
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Appendix A 

Derivation of Coefficients of System Responses 

To find the polynomials of the input current it is simply a matter of employing 
Ohm’s Law in the frequency domain, 


i Ci (s) 


Z. q (s) ' 


When the Mapham inverter is connected to a power system which can be modelled as 
a series R-L circuit, the equivalent input impedance of the entire network is the impedance 
shown in Figure A.1 . The following steps are used to find Z eq which is connected to 
the source. 



Figure A.1 - Total System Connected to the Source. 


1. Series of sL eq) sL eqI and R eq yields Z^s). 

Z 1 ( s) = R eq + S ( Lj + L^g) . 
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2. Parallel of Z.,(s) and sL m yields Z 2 (s). 


Z 2 (s) 


5J_S) ( SLJ 
Z 2 (s) + SL m 


s*L m U,*L„) * sR„L, 
s(Lj 


3. Series of Z 2 (s) and 1 /sC s yields Z 3 (s). 


Zy ( S ) = + Z a ( S ) 

3 SC c * 




sC„ 


siL^ L eq + L m ) + 22, 


e<y 


_ +s 2 i? ag L„C s + s(L J + L. q +L iB ) + 22, 

S 2 C a ( L x + L gq + L m ) +SR. 9 C e 


4. Parallel of Z 3 (s) and 1 /C r yields Z 4 (s). 


z 4 (s) = 


.tel 


( s ) 


hr + z 3< s > 


SC, 


S*C s L m (L 1 + L e<1 ) +S 2 R eq L m C s +sU 1 + L il _ + L n ) + 22 


_£ 2 _ 


S* C s C l L m (L 1 +L eq ) + S 3 L. Cj, C s + S 2 C g ( Lj + + L m ) + S22 e(? (C z + C a ) 
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5. Series Z 4 (s) and sL r yields Z eq (s). 


Z^(8) « SL Z * Z a (8) 


j»L,t.C < C J <ti*Z. w ,) ♦Z. J (C,«C r ) Uj »£„♦!,» 1 ♦ »»<«,. U, <<?,♦ C r ) ♦ L.C,H ♦ «<L, ♦ 

**L J1 C.e,<L J ♦£„) [<C,»C,> <t|*L^.l.>) ♦»*„<C,»C,> 


Since the Laplace transform of the switched input voltage is 




then the input current is 


I e,<«) 



(s) 


s J L m C $ C t iLj + L^) ^S 1 L m C t C t R^*s[ lC t + C M ) (Ij » I w * E«) 1 »it #q (C < »C f ) 


# S I,,Z,.C # C X Uj + L.*) ♦ ff* [ L.C, ( Lj ♦ L^) f ♦ X..) ] ♦ M t < C 9 + C x ) ♦ I,C,1 > ♦ # (Ij ♦ I.) 


The denominator of the above equation yields the coefficients which are necessary to find 
the natural response of the input current. 

When the power system equivalent is a series R-C circuit rather than the series 
R-L which is shown above, the method outlined above is used to find the coefficients of 
the polynomial in the denominator of the current's function in the frequency domain. 
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